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We present the first QCD spectral sum rules analysis of the SU (3) breaking parameter £ and an improved estimate 
of the renormalization group invariant (RGI) bag constant Bs q both entering into the S -B^ s mass-differences. 

The averages of the results from the Laplace and moment sum rules to order a s are fg \J Bb — (247 ± 59) MeV 

and £ = fs a y/ Bb s //s \pB~B — (1.18 ± 0.03), in units where = 130.7 MeV. Combined with the experimental data 
on the mass-differences AM^,, one obtains the constraint on the CKM weak mixing angle \V ts /Vtd\ 2 > 20.0(1.1). 
Alternatively, using the weak mixing angle from the analysis of the unitarity triangle and the data on AMj, one 
predicts AA1 S = 18.6(2.2) ps _1 in agreement with the present experimental lower bound and within the reach of 
Tevatron 2. 
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1. Introduction 

B? s and B? , are not eigenstates of the weak 

(s) (s) t> 

Hamiltonian, such that their oscillation frequency 
is governed by their mass-difference AMq. The 
measurement by the UA1 collaboration |jj of a 
large value of AMd was the first indication of the 
heavy top quark mass. In the SM, the mass- 
difference is approximately given by |2| 



AM„ 



^M^\V tq V t * b \ 



So 
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(B q \O q (x)\B° q ) 
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where the AB = 2 local operator O q (x) is defined 
as 



O q (x) = (b^Lq)(b^Lq) , 



(2) 



with L = (1 — 75)/2 and q = d,s. So, f]B a- n d 
Cb (y) are short distance quantities calculable per- 
turbatively. (Here we are following the notation of 
Ref. S.) On the other hand, the matrix element 
(B q \O q \B^) requires non-perturbative QCD calcu- 



lations, and is usually parametrized for as 



(B° q \O q \B° q ) 
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/b is the B q decay constant normalized as f„ = 
130.7 MeV, and Bs q is the so-called bag parameter 
which is Bs q = 1 if one uses a vacuum saturation 
of the matrix element and equal to 3/4 in the large 
N c limit. From Eq. ([!]), it is clear that the mea- 
surement of AMd provides one of the CKM mixing 
angles \Vtd\ if one uses |Vt&| ~ 1. One can also 
extract this quantity from the ratio 
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(4) 



since in the SM with three generations and unitar- 
ity constraints, |V* S | — \V c b\- Here 



(5) 



The great advantage of Eq. (|J) compared with 
the former relation in Eq. (Ill) is that in the ra- 
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tio, different systematics in the evaluation of the 
matrix element tends to cancel out, thus provid- 
ing a more accurate prediction. However, unlike 
AMd — 0.479(12) ps^ 1 , which is measured with 
a good precision || , the determination of AM S is 
an experimental challenge due to the rapid oscilla- 
tion of the B®-B® system. At present, only a lower 
bound of 13.1 ps _1 is available at the 95% CL from 
experiments || , but this bound already provides a 
strong constraint on \ Vtd\- 

2. Two-point function sum rule 

Ref. Q has extended the analysis of the K°-K° 
systems of Ref. M, using two-point correlator of 
the four-quark operators into the analysis of the 
quantity f B \/B B which governs the B°-B° mass 
difference. The two-point correlator defined as 

^ H { q 2 ) = ijd 4 x (0\TO q (x) (0,(0))* |0> , (6) 

is built from the AB — 2 weak operator given 
in Eq. (Q). The two-point function approach is 
very convenient due to its simple analytic proper- 
ties which are not the case of approach based on 
three-point functions Q However, it involves non- 
trivial QCD calculations which become technically 
complicated when one includes the contributions 
of radiative corrections due to non-factorizable di- 
agrams. These perturbative (PT) radiative correc- 
tions due to factorizable and non-factorizable dia- 
grams have been already computed in Ref. (re- 
ferred as NP), where it has been found that the 
factorizable corrections are large while the non- 
factorizable ones are negligibly small. NP analysis 
has confirmed the estimate in Ref. from lowest 
order calculations, where under some assumptions 
on the contributions of higher mass resonances to 
the spectral function, the value of the bag param- 
eter Bb has been found to be 



B Bd (4M b 2 )^(l±0.15) 



(7) 



This value is comparable with the one B Bd = 1 
from the vacuum saturation estimate, which is ex- 
pected to be a quite good approximation due to the 
relative high-scale of the B-meson mass. Equiva- 
lently, the corresponding RGI quantity is 



B Bd ~ (1.5 ±0.2), 

where we have used the relation 

B B =B B (u)a?!l + (^-)(^ 

1 For detailed criticisms, sec Iql. 



(8) 



(9) 



with 70 = 1 being the anomalous dimension of the 
operator O q and /3\ — —23/6 for 5 flavours, v is 
the subtraction point. The NLO corrections have 
been obtained in the MS scheme § . We have also 
used the value of the bottom quark pole mass fipfl 



M b = (4.66 ± 0.06) GeV 



(10) 



In the following, we study (for the first time), from 
the QCD spectral sum rules (QSSR) method |, the 
SU (3) breaking effects on the ratio £ defined previ- 
ously in Eq. (^|), where a similar analysis of the ra- 
tios of the decay constants has given the values |1Q| 



1.15 ±0.04, ^ 

JD }b 



1.16 ±0.04. 



(11) 



We shall also improve the previous result of Ref. [Q, 
|?j on B Bd by the inclusion of the B q B* and B*B* 
resonances into the spectral function. 

3. Inputs for the sum rule analysis 

We shall be concerned here with the two-point cor- 
relator defined in Eq. (^) . The hadronic part of the 
spectral function can be conveniently parametrized 
using the effective realization M 



1 



Of=-g g d, 1 B q d»B q 



(12) 



where 



f B , B B a 



indicates higher resonances and g q = 
i B ±j Bq . Retaining the BB, BB* and B*B* res- 
onances and parametrizing the higher resonances 
with the QCD continuum contribution, it gives 



-lmiPf(t) 

7T 




0{t - (M 



(13) 



below the QCD continuum threshold t c . The func- 
tion X(x, y, z) is a phase space factor, 

A(a;, y, z) = x 2 + y 2 + z 2 - 2xy - 2yz - 2zx. (14) 

2 Our preliminary results have been presented in . 
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We have used, to a first approximation, the large 
Mb and vacuum saturation relations: 



9q = fl3 q B B q ^ 9q' = /bj 



(15) 



among the couplings. The results fs ~ /s* have 
been also obtained from QCD spectral sum rules 
(|], while the vacuum saturation relation Bb* — 
1 ~ Bs q is a posteriori expected to be a good 
approximation as indicated by the result obtained 
later on in this paper. The short distance ex- 
pression of the spectral function is obtained using 
the Operator Product Expansion (OPE) including 
non-perturbative condensates Q 0. The massless 
(m q — 0) expression for the lowest perturbative 
and gluon condensate contributions has been ob- 
tained in Ref. Q . Radiative factorizable and non- 
factorizable corrections to the perturbative graphs 
in the massless light quark case have been obtained 
in NP §. 

4. SU(3) breaking contributions 

The lowest order perturbative contribution for 
m s 7^ to the two-point correlator is 

-K erl (i) - 6(t - 4(M 6 + m s f) 

7T 
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Here 8 = Mjf/t and 8' = m 2 /t, respectively. The 
functions f(x,y) and g(x,y) are defined by 

f{x,y) = 2-x-y-(x-y) 2 , (17) 



j(x, y) = 1 + x + y - 2{x - y f 



(18) 



We include the 0(a s ) correction from factorizable 
diagrams by using the results in the MS scheme 



3 We shall not include here the effects of tachyonic gluon 
mass [ha] or some other renormalon-like terms which give 



small effects in various examples |13l|l 



for the two-point correlators of currents [jl5j [] This 
can be done using the convolution formula, 

-ImC* (t) = 9(t - 4(M 6 + m s f ) 

IT 

X S"T / _ dz / _ du 

07r JiVs+Vs 7 ) 2 JiVs+Vs 7 ) 2 

xA 1/2 (l,z,u) 

x {lmn^(zi)Imn Qs ^(trf) 

+ lmU^ v (zt)lmU 0flv {ut)} . (19) 

Here U^^q 2 ) and U^(q 2 ) are respectively the low- 
est and the next-to-leading order QCD contribu- 
tion to the two point correlator (q 2 ) defined by 



x{0\T(b(x)^Ls(x))(s(0)^Lb(0))\0). 

The quark condensate contribution reads 
1. 



(20) 
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IT 

^-4(M b + m s ) 2 )- 384 ^... s _ / 

X dq 2 1 ^U + 2q 2 — 

J(M b +m s ) 2 V Oq 

q q 



M 2 



$)l<f-.v£->iV 



Here Aq, Ai, and /i are defined by 



Ao = A ( 1, — 

q z q A 



Ai = A ( 1, — j-, —4- ) , 



Ti 11 



Mr 



2 >2 



h = 1 + ^ + ^-2 



41 



41 



3i 



■(21) 

(22) 
(23) 
(24) 



5. The sum rule analysis 

For the sum rule analysis, we shall work like 
with the moments 



M { q n) = / dt t n -lm%{t) , 



(25) 



4 We shall neglect the nonfactorizable corrections in our anal- 
ysis, according to the results in NP |7| obtained in a slight 
variant of the MS scheme. 
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Figure 1. Moment sum rules analysis of /S( 3 ) y B (a) 

for different values of r q and n: /bV versus: a) r d at 

n = —30, b) n at r d = 1.13; fs E \/ Bs a versus: c) r s at 
n = —26, d) versus n at r s = 1.17. Dotted curve: lowest 
order perturbative contribution; dashed curve: lowest or- 
der perturbative + m s (ss) [only for c) and d)] + (ct s G 2 ) 
condensates. Solid curves: total contribution to order ct s . 



and with the Laplace sum rule 

, to) 

1. 



dt 



4(M b +m q ) 



e " —Imip q (t) 

7T 



(26) 



In so doing, in addition to the pQCD input param- 
eters given previously, we shall need the values of 
the QCD condensates and SU(S) breaking param- 
eters, which we give in Table 1. We show in Fig. [I] 

the moment sum rules analysis of /b (s) \J Ab (s) for 
different values of r q and n. As one can see from 
Fig. [l]a, b, the stability regions of the quantity, 



>B,, 



(27) 



versus the number n of moments and the contin- 
uum threshold 

Ad) 

rd = ^7o , (28) 



4M fe 2 ' 

are obtained for large ranges ending with an ex- 
tremum for 



-30 



1.13 . 



(29) 



These range of values of the sum rule parame- 
ters are in good agreement with previous results 
in Ref. Q and NP [Q. Analogous values of n and 



r s stabilities are also obtained in the analysis of 

VFs = Ib 3 \J B Bb (see Fig. [l]c, d), with 

n~-26, r s ~ 1.17 . (30) 

One can notice that the stabilities in the contin- 
uum for gd and g s differ slightly as a reflection of 
the SU(3) breakings, which one can parametrize 
numerically as 

r s ~ «r d + 0.05 . (31) 

Similar analysis is done with the Laplace sum rules. 



We show in Fig. || the predictions of }b {b) y Bb {s) 
for different values of r q and r, where an extremum 
is obtained for 



r ~ 0.3 GeV" 



(32) 
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Figure 2. Laplace sum rules analysis of fB, 3 s \ Bb^ 



for 



different values of r q and r: /bV Bg versus: a) r,j at r = 
0.31GeV -2 , b) t at r d = 1.13; f Ba y/B Bs versus: c) r s at 
r = 0.26GeV~ 2 , d)versus r at r s = 1.20. The curves are 
the same as in Fig. Ill 



6. Results and implications on 
and AM, 



\V ta /Vt 



td\ 



We take as a conservative result for gd, from the 
moments sum rule analysis, the one from a large 
range of n = —10 to —30 and for rd — 1.06 to 1.17. 
Adding quadratically the different sources of errors 
in Table 1, we obtain 



(245 ± 57) MeV , 



(33) 
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Table 1 

Different sources of errors in the estimate of fg ■> 
that the error is zero or negligible. 



and £ in units where f w = 130.7 MeV. The box marked with — means 



Sources 



[MeV] 



A£ x 10 2 



-n ~ (30 - 10) 
r ~ (0.1 - 0.31) GeV" 2 
r d ~ 1.06 - 1.17 
A 5 = (216±H) MeV 1,@ 
v = M b ~ 2Af(, 
a 2 : geometric PT series 
M 6 = (4.66 ± 0.06) GeV § 
(a s G 2 ) = (0.07 ±0.01) GeV 4 
(mi) (2) = -(254 ± 15)3 MeV 3 
(ss)/(uu) =0.7 ±0.2 [_" 
m s (2) = (117 ±23) MeV 




Moments 


Laplace 


Moments 


8.3 




1.5 




12.4 




7.9 


7.8 


1.0 


0.4 


0.4 


0.1 


8.7 


9.1 


0.2 


43.0 


45.0 


0.6 


34.6 


37.1 


0.9 


1.3 


1.6 








0.2 






0.3 






1.5 



Laplace 



1.7 
2.5 
0.1 
0.2 
0.7 
1.2 

0.3 
0.4 
1.7 



Total 



57.1 



60.8 



2.6 



in units where /„■ = 130.7 MeV. The most rele- 
vant errors given in Eq. ( |33"|) come from Mb and 
the truncation of the PT series. We have es- 
timated the latter by assuming that the coeffi- 
cient of the a 2 contribution comes from a geomet- 
ric growth of the PT coefficients. The other pa- 
rameters n, r"d,A, (qq), (a s G 2 ) and v (subtraction 
point) induce smaller errors as given in Table 1. 
We proceed in a similar way for g s . Then, we take 
the range n = —10 to —26 and r s — 1.10 to 1.21, 
and deduce the ratio 



Then, we deduce 



1.174 ±0.026 



(34) 



JbVBb 

where the errors come almost equally from 
n, r q , m s , Mb and the a 2 term. As expected, 
we have smaller errors for the ratio £ due to the 
cancellation of the systematics. 
We proceed in the same way with the Laplace sum 
rules where we take the range of t values from 0.1 
to 0.37 GeV -2 (see Fig. || a, b) in order to have a 
conservative result. Then, we deduce 




(249 ±61) MeV , 
1.187 ±0.038 . 



(35) 



As a final result, we take the arithmetic average 
from the moments and Laplace sum rules results. 




(247 ± 59) MeV, 
1.18 ±0.03 , 



(36) 



in the unit where f„ = 130.7 MeV. These results 
can be compared with different lattice fg\J Bg ~ 
(230±32) MeV, £ ~ 1.14±0.06, and global-fit of the 

CKM mixing angles giving fg\J Bg — (231 ± 15) 
MeV quoted in |l^ , po|| . By comparing our results 
Eq. (H) with the one of f B JfB d in Eq. @ [J one 
can conclude (to a good approximation) that 



(1.65 ±0.38) 
(1.1 ±0.25) , 



(37) 



indicating a negligible SU (3) breaking for the bag 
parameter. For a consistency, we have used the 
estimate to order a, [E2| 



f B ^ (1.47 ± 0.10)/^ ~ (192 ± 19) MeV , 



(38) 



and we have assumed that the error from fg com- 
pensates the one in Eq. (|3^) . The result is in excel- 
lent agreement with the previous result of Ref. § 

5 One can notice that similar strengths of the SU(3) break- 
ings have been obtained for the B — » K*^y and B — > Klv 
form factors [EM. 



5 



in Eqs. (|7|) and (ph, and agrees within the errors 
with the lattice estimates 
Using the experimental values 



AM d = 0.479(12) ps- 1 , 
AM S > 13.1 ps- 1 (95% CL), 

one can deduce from Eq. (Q) 

2 



> 20.0(1.1) 



Alternatively, using 
1 



Psd 



A 2 [(1 - p) 2 + fj 2 



28.4(2.9) 



(39) 



(40) 



(41) 



with Jl9| 
A ~ 

/' ( 1 
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0.2237(33) , 
A^ 
2 

A 2 
2 



r, 1 



0.223(38) , 
0.316(40) , 



(42) 



A, p and rj being the Wolfenstein parameters, we 
deduce 



AM S ~ 18.6(2.2) ps- 1 



(43) 



in agreement with the present experimental lower 
bound and within the reach of Tevatron run 2 ex- 
periments. 

7. Conclusions 

We have applied QCD spectral sum rules for ex- 
tracting (for the first time) the SU(3) breaking pa- 
rameter £, and for improving the estimate of the 
quantity /b\/ Bb- Our predictions are given in 
Eq. (|36|). The phenomenological consequences of 
our results for the CKM mixing angle and B^ - 
B® mass-differences are given in Eqs. fl4C| ) and 
© • 
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